Isospin-dependent pion in-medium effects on charged pion ratio in heavy ion collisions 
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Using results from the chiral perturbation theory for the s-wave interaction and the delta- 
resonance model for the p-wave interaction of pions with nucleons, we have evaluated the spectral 
functions of pions in asymmetric nuclear matter with unequal proton and neutron densities. We 
find that in hot dense neutron-rich matter the strength of the spectral function of positively charged 
pion at low energies is somewhat larger than that of negatively charged pion. In a thermal model, 
this isospin-dependent effect slightly reduces the ratio of negatively charged to positively charged 
pions that are produced in heavy ion collisions induced by radioactive beams. Relevance of our 
results to the determination of the nuclear symmetry energy from measured ratio of negatively to 
positively charged pions produced in heavy ion collisions is discussed. 

PACS numbers: 21.65.Jk, 21.65.Ef, 25.75.Dw 
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I. INTRODUCTION 

The nuclear symmetry energy is the energy needed per 
nucleon to convert all protons in a symmetric nuclear 
matter to neutrons. Knowledge on the density depen- 
dence of the nuclear symmetry energy is important for 
understanding the dynamics of heavy ion collisions in- 
duced by radioactive beams, the structure of exotic nuclei 
with large neutron or proton excess, and many important 
issues in nuclear astrophysics [H-Ilj. At normal nuclear 
matter density, the nuclear symmetry energy has long 
been known to have a value of about 30 MeV from fit- 
ting the binding energies of atomic nuclei with the liquid- 
drop mass formula. Somewhat stringent constraints on 
the nuclear symmetry energy below the normal nuclear 
density have also been obtained during past few years 
from studies of the isospin diffusion and isoscal- 

ing Q in heavy-ion reactions, the size of neutron skin 
in heavy nuclei [Toj], and the isotope dependence of the 
giant monopole resonances in even- A Sn isotopes [TTj | . 
For nuclear symmetry energy at high densities, transport 
model studies have shown that the ratio of negatively to 
positively charged pions produced in heavy ion collisions 
with neutron-rich nuclei is sensitive to its stiffness 
EH, and there has been an attempt [l5[ to extract the 
symmetry energy at supra-saturation densities from the 
FOPI data on the ir~ /ir + ratio in central heavy-ion col- 
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lisions at SIS/GSI [Ig. 

The transport model used in Ref. [H, EH, EH neglects, 
however, medium effects on pions, although it includes 
those on nucleons and produced A resonances through 
their isospin-dependent mean-field potentials and scat- 
tering cross sections. It is well-known that pions interact 
strongly in nuclear medium as a result of their p-wave 
interactions through couplings to the nucleon-particle- 
nucleon-hole and delta-particle-nucleon-holc (A-holc) 
excitations, leading to the softening of their dispersion 
relations or increased strength of their spectral functions 
at low energies |17h22| . Furthermore, it has been shown 
in Ref. (23j that in asymmetric nuclear matter with dif- 
ferent proton and neutron fractions, these medium effects 
depend on the charge of the pion. Including the p-wave 
pion medium effects in the transport model has previ- 
ously been shown to enhance the production of low en- 
ergy pions in high energy heavy ion collisions, although 
it does not affect the total pion yield [24], HH . Properties 
of pions in nuclear medium are also modified by their 
s-wave interactions with nucleons. Studies based on chi- 
ral perturbation theory have shown that in asymmetric 
nuclear matter this effect also depends on the charge of 
the pion (26|, leading to different in-medium masses for 
negatively and positively charged pions. Since pions of 
different charges are modified differently in asymmetric 
nuclear matter, it is of interest to study how such isospin- 
dependent medium effects would affect the ratio of nega- 
tively to positively charged pions produced in heavy ion 
collisions. 

This paper is organized as follows. We first evalu- 
ate in Sec. [TT] the spectral functions of pions and delta- 
resonances in hot dense asymmetric nuclear matter. Us- 
ing these spectral functions, we then study in Sec. lIIII the 
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charged pion ratio in a thermal model and compare the 
results with available experimental data. We conclude 
with a summary in Sec. Qv] 



II. MEDIUM EFFECTS 

A. Pion in-medium properties 

Pions in nuclear medium acquire self-energies through 
their s— wave and p-wave interactions. In asymmetric 
nuclear matter, the pion self energy further depends on 
the charge of the pion. 



1. Pion-nucleon s-wave interaction 

The contribution of the pion-nucleon s-wave inter- 
action to the pion self-energy has been calculated in 
Ref. |26[ up to the two-loop order in chiral perturbation 
theory. In asymmetric nuclear matter of proton density 
p p and neutron density p n , the resulting pion self-energy 
depends on the charge of the pion, and for n~ , 7r + , and 
7T°, their self-energies due to the s-wave interactions are 
given, respectively, by 

^(pp, Pn) = Pn[T~ N - T+ N ] - p p [T- N + T+ N ] 

+ I Vcl(/ 3 f» Pn) + n ™rOp: Pn) 
n+(ft,,Pn) = Uj(Pn,Pp) 

tt°s(Pp,Pn) = ~{Pp+ Pn)T+ N +Il cm {pp,p n ). (1) 

In the above, are the isospin-even and isospin-odd 
■kN s-wave scattering T-matrices, which are given by the 
one-loop contribution in chiral perturbation theory and 
have the empirical values w 1.847 fm and T^. 



nuclear medium at temperature T, its self-energy due to 
the p-wave interaction is then given by 

IC « U^-Yk 2 F*(k) ]T \(lm T \lm t \m T )\ 2 

r d 3 P 1 

X J (27T) 3 e ( m iv+P 2 / 2m N + U N T -" B -2m TfJ , Q )/T ~ 1 
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In the above, m A ~ 1232 MeV is the mass of A res- 
onance; / A ~ 3.5 is the wNA coupling constant and 

i^(Jfc) = [1 + 0.6(fc 2 /™^)]" 1/2 HI is the nNA form fac- 
tor determined by fitting the decay width ~ 118 MeV of 
A resonance in free space. The summation in Eq. is 
over the nucleon isospin state m T and the A resonance 



isospin state m T \ and the factor 



|1 m t \ m T ) is the 



Clebsch-Gordan coefficient from the isospin coupling of 
pion with nucleon and A resonance. The momentum 
integration is over that of nucleons in the nuclear mat- 
ter given by a Fermi-Dirac distribution with p B and pq 
being, respectively, the baryon and charge chemical po- 
tentials determined by charge and baryon number con- 
servations; U 7 ^ is the mean-field potential of nucleons of 
isospin state m T in nuclear matter; and T A T and U* T 
are, respectively, the width and mean-field potential of 
A resonance of isospin state m T . 

For the nucleon mean-field potential U^ T , we use that 
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-0.045 fm extracted from the energy shift and width of from the moment urn- independent (MID) interaction [1|, 



the Is level in pionic hydrogen atom. The term IT^j is 
due to the relativistic correction, whereas the terms n,r ol ., 
and IlJ? or are the contributions from the two-loop order 
in chiral perturbation theory. Their expressions can be 
found in Ref. [26| , where it is also shown that in nuclear 
matter of density p = 0.165 fm~ 3 and isospin asymmetry 
^ = (Pn ~ Pp)/(Pn + Pp) = 0.2, such as in the interior 
of a Pb nucleus, changes of the pion masses due to their 
s-wave interactions are Am,- = 13.8 MeV, Am 7r + = 
-1.2 MeV, and Am,» = 6.1 MeV. Compared to their 
masses in free space, the tt~ mass increases whereas the 
7r + mass decreases in neutron-rich nuclear matter. 



2. Pion-nucleon p-wave interaction 

For the p-wave interactions of pions in nuclear matter, 
we consider only the dominant A-hole excitations as in 
Ref. [13], as the contribution from the nucleon particle- 
hole excitations is known to be small. For a pion of 
isospin state mt, energy oj, and momentum k in a hot 
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being the nucleon symmetry potential. The parameters 
a = -293.4 MeV, (3 = 240.1 MeV, and 7 = 1.216 
are chosen to give a compressibility of 212 MeV and a 
binding energy per nucleon of —16 MeV for symmet- 
ric nuclear matter at saturation or normal nuclear den- 
sity p Q = 0.16 fm -3 . The nucleon symmetry poten- 
tial U™j (p B , <5iikc) depends on the baryon density p B = 
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P n + P P + Pa- + Pa" + Pa+ + *°A++ an d tne isospin asym- 
metry <5 likc = {p n -p p + p A --p A+ ++p A o/3-p A+ /3)/p B 
of asymmetric hadronic matter, which is a generalization 
of the isospin asymmetry S usually defined for asymmet- 
ric nuclear matter without A resonances fl2| . The nu- 
cleon mean-field potential also depends on the stiffness 
of nuclear symmetry energy through the parameter x via 
the functions F(x) and G(x). We consider the three 
cases of x = 0, x = 0.5, and x = 1 with correspond- 
ing values F(x = 0) = 129.98 and G{x = 0) = 1.059, 
F(x = 0.5) = 85.54 and G{x = 0.5) = 1.212, and 
F(x = 1) = 107.23 and G(x = 1) = 1.246. The resulting 
nuclear symmetry energy becomes increasingly softer as 
the value of x increases, with x = 1 giving a nuclear sym- 
metry energy that becomes negative at about three times 
normal nuclear matter density as shown in Fig. [1] These 
symmetry energies reflect the uncertainties in the theo- 
retical predictions on the stiffness of nuclear symmetry 
energy at high densities. 
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3. pion spectral function 



and n™ T due 



In terms of the pion self energies IT 
to the pion-nucleon s- and p-wave interactions in nuclear 
medium, the pion in-medium propagator is given by 

'( W > fc ) = ~ p ~ mHi -nmt 7. . j„\ • ( 6 ) 
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The pion spectral function S™*(u>,k) is related to 
the imaginary part of the pion in-medium propagator 
through 



S^(uj,k) = lmD mt (oj,k). 

7T 



(7) 



B. Delta resonance in-medium properties 
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FIG. 1: (Color online) Nuclear symmetry energies as func- 
tions of nuclear density from the MID interaction for different 
values of the symmetry energy parameter x. 



For the mean-field potentials of A resonances, their 
isoscalar potentials are assumed to be the same as those 
of nucleons, and their symmetry potentials are taken to 
be the average of those for neutrons and protons with 
the weighting factors depending on the charge state of A 

T A+ + _ rip ttA+ _ 2, 



resonances [28|, i.e., U~. 
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Including the short-range A-hole repulsive interaction 
via the Mig dal parameter g' , which has values 1/3 < 
g' < 0.6 [17H22| |. modifies the pion self-energy due to the 



The modification of the pion properties in nuclear 
medium affects the decay width and mass distribution of 
A resonances. For a A resonance of isospin state m T and 
mass M and at rest in nuclear matter, its decay width is 
given by [jj 
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In the above, the first term in the last line is due to 
the decay of the A resonance to the pion but corrected 
by the contact interaction at the irNA vertex, while the 
second term contains the contribution from its decay to 
the A-hole state without coupling to the pion. The two 
temperature-dependent factors in the momentum inte- 
gral take into account, respectively, the Bose enhance- 
ment for the pion and the Pauli blocking of the nucleon. 
To include possible chemical non-equilibrium effect, a fu- 
gacity parameter is introduced for pions. The pion 
energy uj is determined from energy conservation, i.e., 
M + U A T = u} + m N + k 2 /2m N + U'™ T . The resulting 
mass distribution of A resonances is then given by 



P A (M)=A- 



\M)/2 



(9) 



(M~m A ) 2 +T' A T (M)/4 

where A is a normalization constant to ensure the inte 
gration of Pa(M) over M is one. 
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FIG. 2: (Color online) Spectral functions of pions in asym- 
metric nuclear matter of density 2p and isospin asymmetry 
<5iikc = 0.135 as functions of pion energy for different pion 
momenta of (a) m^, (b) 2m,,-, (c) 3m w , and (d) 4m T . All are 
calculated with the Migdal parameter g' — 1/3. 



FIG. 3: (Color online) Mass distributions of A resonances at 
rest in asymmetric nuclear matter of density 2p and isospin 
asymmetry <5iik = 0.135. The solid line corresponds to that 
in free space. The distributions near the threshold and at the 
peak are enlarged in the insets. 



C. Results 

Equations ©, © and flSJ) are solved self-consistently 
to obtain the pion spectral functions and the mass distri- 
butions of A resonances in asymmetric nuclear matter. 
The results obtained with the Migdal parameter g' = 1/3 
are illustrated in Fig. [5] and Fig. [3] for an asymmetric nu- 
clear matter of isospin asymmetry Sy^ ~ 0.135, twice the 
normal nuclear matter density p B = 2p , temperature 
T ~ 43.6 MeV, and chemical potentials [i B ~ 942 MeV 
and Hq ~ —18.5 MeV, corresponding to those to be used 
in the thermal model and also similar to those reached 
in the IBUU transport model with the nuclear symme- 
try energy x = 1 for central Au+Au collisions at the 
beam energy of 0.4 AGeV (IBf . Shown in Fig. [2] are the 
pion spectral functions as functions of pion energy for 
different values of the pion momentum. It is seen that 
for low pion momenta the spectral function at low ener- 
gies has a larger strength for n + (dashed line) than for 
ir° (solid line), which has a strength larger than that for 
7r~ (dotted line), and this behavior is reversed at high 
pion energies. Figure |3] shows the mass distributions of 
A resonances at rest in asymmetric nuclear matter as 
functions of mass. One sees that it is similar to that 
in free space (solid line) as a result of the cancelation 
between the pion in-medium effects, which enhance the 
strength at low masses, and the Pauli-blocking of the nu- 
cleon from delta decay, which reduces the strength at low 
masses. This is consistent with observed similar energy 
dependence of the photo-proton and photo-nucleus ab- 



sorption cross sections around the A resonance [29] . Fur- 
thermore, whereas the strength around the peak of the 
A resonance mass distribution decreases with increasing 
charge of A resonance due to nonzero isospin asymme- 
try of the nuclear medium, that near the threshold in- 
creases with increasing A resonance charge. We note that 
the pion-nucleon s-wave and p-wave interactions have op- 
posite effects on the in-medium properties of pions and 
delta resonances. While the pion-nucleon s-wave inter- 
action increases the tt~ mass and reduces the tt + mass, 
the pion-nucleon p-wave interaction softens the disper- 
sion relation of ir~ more than that of tt + . As a result, 
including only medium effects due to the pion-nucleon 
p-wave interaction would lead to an opposite result, i.e., 
for low momentum pions the spectral function at low en- 
ergies has a larger strength for n~ than for 7r + , although 
the strength around the peak of the A resonance mass 
distribution still decreases with increasing charge of A 
resonance. 



III. CHARGED PION RATIO IN HOT DENSE 
ASYMMETRIC NUCLEAR MATTER 

To see the isospin-dependent pion in-medium effects on 
the ir~ /n + ratio in heavy ion collisions, we use a thermal 
model which assumes that pions are in thermal equilib- 
rium with nucleons and A resonances (30j . The density 
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of particle species i is then given by 



d 3 k 

(2tt) 3 



duj n 'Si(u,k)- 



1 



-l e (u-B^ B -Q iP , Q )/T ±1 - 



(10) 



In the above, g,, Qi, and 2^ are the degeneracy, 
baryon number, charge, and fugacity of the particle. 
The exponent n i is 2 for pions and 1 for nucleons and 
A resonances. For pions, we use the spectral functions 
5™* (a;, k) calculated above. For the spectral functions of 
A resonances, we neglect their momentum dependence 
and thus replace the integration over energy uj by that 
over mass M. The ui in the Fermi-Dirac distribution for 
A resonances is then simply uj = M + k 2 /2M + U^ T . For 
nucleons, their spectral functions are taken to be delta 
functions if we neglect the imaginary part of their self- 



k 2 /2m 



N ' 



energies, i.e., S^(lu, k) — 8{u — m N 

According to studies based on transport models [12| . 
Ell Hil ; the total number of pions and A resonances in 
heavy ion collisions reaches a maximum value when the 
colliding matter achieves the maximum density, and re- 
mains essentially constant during the expansion of the 
matter. For Au+Au collisions at the beam energy of 
0.4 AGeV, for which the tt~/tt + ratio has been mea- 
sured by the FOPI Collaboration at GSI 0, the IBUU 
transport model gives a maximum density that is about 
twice the normal nuclear matter density and is insensi- 
tive to the stiffness of the nuclear symmetry energy, as 
it is mainly determined by the isoscalar part of the nu- 
clear equation of state [l5j. We thus use this density 
in the thermal model. The temperature in the ther- 
mal model is determined by fitting the measured pion 
to nucleon ratio, which is about 0.014 including pions 
and nucleons from decays of A resonances [l6j], without 
medium effects and with unity fugacity parameters for 
all particles, and the value is T ~ 43.6 MeV. The as- 
sumption that pions and A resonances are in chemical 
equilibrium is consistent with the short chemical equili- 
bration times estimated from the pion and A resonance 
production rates. The isospin asymmetry of the hadronic 
matter is then taken to be <5i;k c — 0.080, 0.106, and 0.143, 
corresponding to net charge densities of 0.920p , 0.894p 
and 0.857p , for the three symmetry energies given by 
x = 0, 0.5, and 1, respectively, in order to reproduce 
the 1T-/TT+ ratios of 2.20, 2.40, and 2.60 predicted by 
the IBUU transport model of Ref. [HI using correspond- 
ing symmetry energies without pion in-medium effects. 
Since medium effects enhance the pion and A resonance 
densities, to maintain same pion to nucleon ratio as the 
measured one requires fugacity parameters for pions and 
A resonances to be less than one. Also, the pion in- 
medium effects have been shown to affect only slightly 
the pion and the A resonance abundance [24| , indicating 
that both pions and A resonances are out of chemical 
equilibrium with nucleons when medium effects are in- 
cluded, as expected from estimated increasing pion and 
A resonance chemical equilibration times as a result of 
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FIG. 4: (Color online) The tt~ /ir + ratio in Au+Au collisions 
at the beam energy of 0.4 AGeV for different values of nu- 
clear symmetry energy parameter (x — 0, 0.5, and 1) and 
the Migdal parameter g = 1/3, 0.4, 0.5, and 0.6. Results for 
g' — oo correspond to the case without the pion in-medium 
effects. 



medium effects. Because of the small number of pions 
(about 0.3%) and A (about 1.1%) resonances in the mat- 
ter, the density, temperature, and net charge density of 
the hadronic matter are expected to remain unchanged 
when pion in-medium effects are introduced. They lead 
to, however, a slight reduction of the isospin asymmetry 
to <5iiko — 0.073, 0.098, and 0.135 for the three symmetry 
energies, given by a; = 0, 0.5 and 1, respectively. We note 
that with the fugacity of nucleons kept at zjv = 1, the fu- 
gacity parameters of about z„ = 0.0855 and z A = 0.459 
are needed for the symmetry energy given by x = 1 to 
keep the total number of pions as well as that of A reso- 
nances the same as in the case without pion in-medium 
effects, and that the required values for the fugacity pa- 
rameters are similar for the other two symmetry energies 
considered here. 

Results on the ir~ /tt + ratio in Au+Au collisions at 
the beam energy of 0.4 AGeV are shown in Fig. 2J With 
the value g' — 1/3 for the Migdal parameter, values for 
the 7T~ /ir + ratio are 2.12, 2.30, and 2.47 for the three 
symmetry energy parameters x = 0, 0.5, and 1, respec- 
tively, which are slightly smaller than corresponding val- 
ues for the case without including the pion in-medium 
effects as shown by those for g' ~ oo in Fig. Q] The 
measured ir~ ratio of about 3 by the FOPI Col- 
laboration, shown in Fig. 0] by the dash-dotted line to- 
gether with the error bar, which without the inclusion 
of the pion in-medium effects in the IBUU model favors 
a nuclear symmetry energy softer than the one given by 
x = 1 , thus requires an even softer one after including the 
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isospin-dependent pion in-medium effects. Figure [4] fur- 
ther shows results obtained with larger values of g 1 = 0.4, 
0.5 and 0.6 for the Migdal parameter. It is seen that the 
isospin-dependent pion in-medium effects in these cases 
are not much different from the case of g' — 1/3, indi- 
cating that the pion abundance in hot dense asymmetric 
nuclear matter is affected more by the s-wave, which does 
not depend on g' ', than by the p-wave interaction between 
pion and nucleon. 

IV. SUMMARY 

In summary, we have studied the dependence of the 
pion spectral function in asymmetric nuclear matter on 
the charge of the pion by using results from the chiral per- 
turbation theory for the pion-nucleon s-wave interaction 
and from the A-hole model for the pion-nucleon p-wave 
interaction. Because of increasing n~ and decreasing 7r + 
in-medium masses due to the pion-nucleon s-wave inter- 
action in neutron-rich matter, the strength of tt + spectral 
function at low energies is somewhat larger than that of 
ir~ spectral function, and the strength around the peak 
of the A resonance mass distribution decreases while that 
near the threshold increases with increasing charge of the 
A resonance. In a thermal model that assumes that nu- 
cleons, pions, and A resonances produced in heavy ion 
collisions are in thermal but not chemical equilibrium, 
with the latter needed to maintain the final pion to nu- 
cleon ratio, the ir~ /ir + ratio is slightly reduced in com- 
parison with the case without pion in-medium effects. 
Taking into consideration of the isospin-dependent pion 



in-medium effects in the transport model thus would have 
some, albeit not very significant, effects on the extraction 
of the nuclear symmetry energy from measured tt~ /n + 
ratio. On the other hand, the predicted 7r _ /7r + ratio for 
a given parametrization of the symmetry energy can be 
quite different in different transport models. For exam- 
ple, it was shown in Ref. [l4| that the ir~ /tt + ratio in hot 
dense nuclear matter increases with increasing stiffness 
of the nuclear symmetry energy at high densities if the 
latter is parameterized using the relativistic mean-field 
model, opposite to the results from the IBUU model [H[ 
which shows a decreasing 7r~ /n + ratio with increasing 
stiffness of the nuclear symmetry energy. Also, it was re- 
cently shown in an improved isospin-dependent quantum 
molecular dynamic (ImlQMD) model that the measured 
7r~ /ir + ratio from the FOPI collaboration is consistent 
with a nuclear symmetry energy that is even stiffer than 
the one corresponding to the symmetry energy parame- 
ter x — in the MDI interaction [3 Incontradictory to 
the conclusion from the IBUU model [15j that a symme- 
try energy softer than x = 1 is needed to reproduce the 
measured ?t~ /n + ratio. Further theoretical work is thus 
needed to understand the relation between the tt~ /n + 
ratio and the behavior of the nuclear symmetry energy 
at high densities in the transport model description of 
heavy ion collisions. 
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